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Résumé 

Nous considérons un opérateur /i-pseudodifférentiel non- autoadjoint 
dans la limite semiclassique. p désigne le symbole principal. Nous sa- 
vons que la résolvante existe à l'intérieur de l'image de p jusqu'à une 
distance 0({h\n ji) k+1 ), de certains points du bord, où k £ {2,4, . . .}. 
Dans ce travail, nous précisons les estimations de résolvantes qu'ont 
obtenues différents auteurs dans le cas k = 2, et en dimension 1. Pour 
la preuve, il s'agit de construire, via un scaling, des quasimodes pour 
des valeurs du paramètre spectral très proches du bord de l'image de 
p. 

Abstract 

We consider a non-self-adjoint pseudodifferential operator in the 
semi-classical limit (h — > 0). The principal symbol is given by p. 
We know that the résolvent {z — P) -1 exists inside the range up to 

a distance 0((/iln ¥) k + 1 ) from certain boundary points, where k E 
{2,4, . . .}. In this work, we improve the résolvent estimâtes given by 
différent authors in the case k = 2, and in dimension one. For the 
proof, we will construct quasimodes by a scaling for z very close to 
the boundary. 
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1 Introduction 

Dans [5], N. Dencker, J. Sjostrand et M. Zworski obtiennent une estima- 
tion de la résolvante pour certains points du bord de l'image du symbole 
principal. Dans le cas d'un point zq de type fini d'ordre k, ils obtiennent 

k 

que la résolvante a une croissance en h fc +! , dans des disques de rayon 
0(/i*+ï) centrés en z . Dans [16], J. Sjostrand montre que la résolvante peut 
s'étendre à des disques de rayon 0((h ln |) *+ î ), et donne une majoration 
pour la résolvante. Une majoration similaire avait été obtenue auparavant 
par l'auteur pour un opérateur modèle hD x + g(x), x G S 1 , avec k = 2, voir 
PQ Chapitre 4. On citera aussi J. Martinet [12] qui obtient un encadrement 
de la norme de la résolvante pour l'opérateur d'Airy complexe. La preuve 
est basée sur une analyse directe du semigroupe après conjugaison par la 
transformée de Fourier. Dans ce travail, il s'agit d'améliorer et généraliser les 
estimations de la résolvante obtenues par ces différents auteurs. On se limite 
au cas d'un point d'ordre 2, et à la dimension 1. Les estimations obtenues 
nous permettent de préciser l'estimation de résolvante obtenue par Martinet 
[T2"] pour l'opérateur d'Airy complexe sur la droite réelle (D x ) 2 + ix, i G 1, et 
d'obtenir le comportement à l'infini des lignes de niveaux de la résolvante de 
l'oscillateur harmonique non- autoadjoint. Ce dernier résultat vient compléter 
ceux publiés par L. Boulton [2], K. Pravda-Starov |15] . Mentionnons aussi 
le travail [3] de Davies et Kuijlaars qui obtiennent le premier terme dans 
l'asymptotique du logarithme de la norme des projecteurs spectraux pour les 
grandes valeurs propres pour l'oscillateur harmonique non-autoadjoint. 
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Cas modèle sur H. à paramètre. Nous considérons l'opérateur modèle 
non-autoadjoint dans L 2 (R) 

1 d 

P a = hD x + g(x,a) + h 2 g*(x,a;h), h E (0,1], Dx = ~^ 

g et g* dépendent de manière C°° de x et d'un paramètre a strictement 
positif et borné. On pose g(x, a; h) := g(x, a) + h 2 g*(x, a; h) et on demande 
qu'il existe C > tel que 



(cette hypothèse peut facilement être affaiblie). Pour g* nous avons le dévelop- 
pement asymptotique en puissances de h suivant 



g*(x, a; h) ~ g®(x, a) + hgl(x, a) + . . . dans C™(R x [0, 1]). (1.1) 



L'opérateur P a est muni du domaine H SC (R) := {u G L 2 (R)\ \\x 2 u\\ + 
||/iD x ti|| < +oo}. On note p a le symbole principale semiclassique de P a , 
soit p a = £ + g(x, a). Rappelons que la définition du crochet de Poisson est 



pour deux fonctions a(x,£), b(x,£) de classe C 1 (M 2 ). Ici H a = (a^d x — a x d^) 
désigne le champ Hamiltonien. 

Proposition 1.1 Soient g comme ci-dessus avec g(x,a) G C°°(Rx]0, C[), 
et pour tout h fixé g*(x, a; h) G C°°(IR x [0, 1/C]) pour un C > donné. Nous 
rajoutons les hypothèses suivantes : 
pour x 7^ 0, g satisfait lmg(x, 0) > 0, et 



L'opérateur P a est muni du domaine H SC (R). Dans ces conditions, la résolvante 
est définie pour tout a 



Reg = 0(x 2 + 1), 



{a,b}(x,0 := (4*4 - a' x b'ç)(x,0 = H a b, 



(1.2) 



<?(0,0) = 0, g' a (0,0) = -i, 
g' x (0,0) = 0, lm<&(0,0)>0. 



(1.3) 
(1.4) 



(hD x + g(x, a; h))' 1 : L 2 (R) ->■ L 2 (R), 



(1.5) 



et vérifie pour h <C o?l 2 et a assez petit, 
iin-in v/7rexp(ilm£ (a)) 




(1 + <9(/i)) + 0( 



1 



^ 1/2 (è{^^«}(p+))^(è{p«'^}(p-))^ 



(1.6) 



3 



où h = h/a 3 / 2 , et p±(a) := (x±(a),Ç±(a)) sont les solutions de p a (p±) = 
£± + ç( x ±, a ) — avec ^Im.'(x±,a) > (impliquant Im g (x±, a) = 0). i es t 
donné par 

£ :=-lm ! g{y,a)dy^a 3 ' 2 . (1.7) 

•/ X-L. 



Nous avons pour a assez petit 



i 



1/2 



VIpo, 2î{po,po}|(o,o)y 

Voir la section 2 pour la preuve. 

Remarque 1.2 Le dénominateur du premier de \l.Œ est de V ordre de s/ha 1 / 4 
et ce qui le rend dominant dès que a h 2 / 3 . 

Remarque 1.3 La proposition reste valable pour a G [0, 1] si on rajoute l'hy- 
pothèse (automatiquement vérifiée pour a petit) que l'équation lmg(x, a) = 
admet exactement deux solution x = x+(a), x_(a) et de plus ^fïmg' x (x±(a), a) 
> 0. 

Nous considérons l'opérateur d'Airy complexe sur la droite réelle A = D 2 +ix, 
étudié par J. Martinet [12], et muni du domaine D(A) = {u G iT 2 (K.) | xu G 
L 2 (M)}. Nous savons que le spectre est vide, que la norme de la résolvante 
(A — z)~ x ne dépend que de Rez (voir Helffer [H]), et qu'en conjuguant par 
une transformée de Fourier, on peut se ramener à l'étude de l'opérateur 

D x — ix 2 + iz. 

Pour z dans la zone d'intérêt Re z > 0, le changement de variable x = 
(Rez)~2y, donne l'opérateur (Rez)Q, où Q est donné par 

1 _ . 2 . Imz 
-Jj T — ix + % — 



(Rez)t Rez 
De la proposition au-dessus et la remarque 11.31 avec h = — ^—r et a = 1, on 

(Rezp 

déduit le corollaire suivant : 
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Corollaire 1.4 Pour Rez > assez grand, la résolvante de A satisfait 
+ix- ~y|(Re Z )-i(l + O(^ipj)) exp(l(Rez)§) 

Ceci précise l'estimation qu'a obtenue J. Martinet dans [12] . La preuve sera 
donnée à la fin de sa sous- section 2.2. 

Cas général en dimension 1. Soit p G S(M. 2 ,m) indépendant de h, où 
m désigne une fonction d'ordre au sens que 

3C > 1, N > tels que m(p) < C (p - p) No m{n), Vp, p G R 2 , 



avec (p — p) = yï + |p — p| 2 . 
L'espace de symboles correspondant est 

» 2 im ) = { a g C°°(M 2 ), |<9>(P)I < C a m(p), p G M 2 , a G N 2 }. 



Nous désignons par P = p w son /i-quantifié de Weyl (voir par exemple jl]), 
que nous considérons dans L 2 (R). Nous faisons une hypothèse d'ellipticité à 
l'infini sur P : 

\p(x,0\>m(x,Ç)/C, |(z,£)| >C, 

m(x,Ç) —> oo, (x,Ç) oo. (1.9) 

Pour garantir que le spectre de P n'est pas le plan complexe, nous avons 
besoin de supposer que 

p(T*M 2 ) C. (1.10) 

Dans ces conditions, le spectre de P est discret pour h assez petit, voir par 
exemple [S] ou [TU] . 

Introduisons l'ensemble 

£(p)=p(T*M), (1.11) 

ce dernier est automatiquement fermé à cause de fil .91) . Nous allons faire une 
hypothèse sur un point Zq G <9£(p) : 

p-\z ) = {p }, {p, {p,p}}(po) ^ 0. (1.12) 

Ceci entraine que dp(p) ^ 0, Nous disons aussi p est de type fini d'ordre 2 
au point z . 
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Le champ Hamiltonien H^^^ est un champ de vecteurs réel et tangent 
à l'ensemble {p| ^{p, p}(p) = 0}. Nous avons 

] - Tt,T [3 s i y p(s) := exp(sH à{ptP} )(p ) C {^{p,P} = 0} 

est une courbe orientée. Sa courbe image par p ; 7 = p o p est une courbe C°° 
avec 7 7^ et dont l'image coïncide avec <9£ près de z . De plus S se situe 
à gauche quand on regarde dans la direction de 7, voir [H]. Près de z , pour 
tout complexe G <9E(p), nous avons p~ l (z*) = {p(s*)}. 7 = p o p est aussi 
une courbe orientée. Le vecteur unitaire tangent à <9£ dans un voisinage de 
zq est alors donné par u{s) := 7(3)/ (7(5) |, lequel vérifie 

/ \ faèfa-P})(p( s )) 
u[s) - \{p,Up>P}Kp( s ))\' 

Près de z , S(p) est un ensemble à bord C°°. Il existe un voisinage W de z 
tel que tout z G S(p) fl W peut s'écrire sous la forme 

z = 7(5) + iau(s), a > 0, z = 7(0). (1-13) 



Proposition 1.5 «Sort P un opérateur h-pseudodifférentiel de symbole indépen- 
dant de h. Nous supposons que (p, z ) vérifie les conditions U.9\) . U.10\) . 
M.lty) . Il existe un voisinage W de z tel que pour tout point z de (S(p) fl 
W) \ dT,(p) nous avons 

p- 1 (z) = { P+ (z),p^z)}, p± = (x±,e±), (1-14) 

avec 

0<±^{p,p}( P± )xa 1 / 2 |p + -p_|xa 1 /2. 

(le point p se scinde en deux points p + et p_ lorsque Von s'éloigne du bord 
deZ{p)).) 

Ici on introduit £o(z). Dans le cas où p^(po) 7^ nous avons la factorisation 

P(p) -z = q(x, £, a, s)(£ + g(x, a, s)) (1.15) 

où si a < alors lmg(x) > 0, x > x_ ou x < x + et Im g(x) < 0, X- < x < x + . 
Ici p±(x±, £-t), £± = — (?(x±, a, s). Sous ces conditions, on définit £q(z) comme 

px-(z) 

£ ( z ) ■=- g(x)dx. (1.16) 

Jx+(z) 
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Si Pç(po) — alors p' x (p) 7^ 0. On applique la discussion précédente kp(x, £) = 
p(—£,x), fonction obtenue de p par composition avec une transformation 
canonique. On pose alors 

î (z) = -f g(x,C)dx (1.17) 

où p(p) — z = q(x, £, a, s)(£ + g(x, a, s)). On peut prouver que 

Proposition 1.6 Si p'ç(po) ^ et p' x (p ) ^ 0, alors nous avons lm£ (z) = 
lm£ (z) + O(a°°). 

Théorème 1.7 II existe une constante T* (< T ,Ti) telle que pour toutes 
constantes C , Ci > il existe une constante C 2 > telle que la résolvante 
(P — z)' 1 est bien définie pour 

h 2 / 3 1 1 

\s\<T„ —^a^C^hln-f 3 , h<—, (1.18) 

et satisfait V estimation 

\\{P - *ni ~ ■ ^^l 11 '" r x (1 + 0(h)) + 0(^— ), 

où h = h/ a 3 / 2 , et £ (qui est une intégrale d'action) vérifie 

£o(z) :— / £cte, (7 relie p_ à p + ). (1-19) 
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à Vienne sponsorisé par K. Groechenig et H. G. Feichtinger et par les grants 
sponsorisés l'un par V. Ivrii et l'autre par M. Sigal. Merci enfin à K. Groe- 
chenig et H. G. Feichtinger à Vienne, et V. Ivrii et M. Sigal. à Toronto pour 
leur accueil. 



2 Modèle sur R à paramètre 

Nous rappelons p a (x, := £ + 9{ x i °) désigne le symbole principal semi- 
classique de P a = hD x + g(x, a; h). Notre opérateur hD x + g(x, a; h) admet 
comme domaine naturel l'espace semiclassique suivant 

H SC (R) := {u G L 2 \ \\u\\h sc '■= \\u\\ + ||^ 2 m|| + ||/i.D x m|| < 00}. 
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L'inverse, définie pour tout a 



{hD x + g(x, a; h))' 1 : L 2 (R) -> L 2 (M) 



(2.1) 



est donnée par la formule suivante 



u(x) 




h 



v(y)dy. 



En particulier, quand a est nul, nous avons (voir [S]) 



\\{hD x + g(x^-h))- l \\ <Ch 



2/3 



(2.2) 



2.1 Ensemble d'énergie pour de petites valeurs du pa- 
ramètre 

Nous cherchons à décrire l'ensemble p^O) pour a <C 1 assez petit. 
D'abord remarquons que 



Résoudre £ + g(x, a) = revient à trouver les solutions de l'équation 
Im g(x, a) = 0. 

Puisque Img^(O) = —1, nous avons par la factorisation lmg(x,a) = 
q(x, a) (a — f(x)), où q ^ dans un voisinage de zéro, avec /(0) = /'(0) = 0, 
et /"(0) = — Img^O) ^ 0. Par le théorème de l'inversion locale, il existe des 
voisinages U + :=]0, &+[ et U- :=]&_, 0[ tels que / soit un C 1 diffémorphisme 
de U + sur ]0, /(&+)[ et de LL sur ]0, /(&_)[. Pour avoir les solutions de 
l'équation a — f(x) = 0, nous inversons simplement une fonction croissante 
en décroissante. Nous trouvons alors 



^(0) = {(0,0)}. 



x G U- 



2 



a = x 2 + 0(\xf) 



(x + 0(x 2 )) 2 . 



Soit 





(2.3) 



Nous déduisons alors le résultat suivant : 



•s 



Lemme 2.1 Soit p a = £ + g(x,a) le symbole principal de P a . U.3\) et 
sont vérifiés. Nous avons, pour a > 0, assez petit, 

Pâ X (0) = {P+(«)} U {/>-(«)}, p±(a) = (x ± ,e±), 

avec 

/ _ 2 \ 1/2 

ï±(a) = T« 5 h ïj - Wnn J + <?(«); (2-4) 

te/ gue 

±^{p a ,^}(p±(a)) >0. (2.5) 

L'ensemble Pg 1 (0) est constitué d'un seul point (0,0) qui va se "scinder" en 
deux p+(a) et p~(a), lorsque a ne s'annule plus. 

Preuve. Il faut remarquer que pour le premier crochet et le second crochet 
de Poisson, nous avons 

Y i iP^P^}(p±) = -Im^(ar±(a),a), 
{Po,7^{Po,Pô}}(0,0) = -lm^(0,0). 

Dès lors, 

±7^{Pa,P^}{p±) = Tlmfl^(a;±(a),a) 

= =F (x±(a)Im <4(0, 0) + almg'^ jO, 0) + 0{x\ + a 2 )) 

(2.6) 

> 0. 

Le premier terme de (12. 6p domine car il se comporte comme y/a. □ 



2.2 Problème de Grushin. 

Nous nous appuyons sur les sections "solutions locales" et "problème de 
Grushin" de Hager [7J et "Enoncé et résolution asymptotique du problème 
de Grushin" de Hager [S] . Pour une lecture "plus aisée" , nous reprenons en 
détail les grandes lignes, et les calculs dont nous aurons besoin. 

On suppose que a > appartient à un intervalle ouvert Q, relativement 
compact, séparé de l'origine par une constante indépendante de h. Nous allons 
d'abord considérer cette situation, puis nous laisserons a devenir petit et nous 
ferons des dilations pour se ramener à ce premier cas. 
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Avec les hypothèses choisies pour g, nous avons 

Va G tt, p-^O) = {p+(a),p-(a)} avec ± j.{p a ,P^}{p±) > 0. (2.7) 
On peut trouver des intervalles J + et J_ disjoints tels que 



x±(Q) (s J±. 

Considérons x± £ C°°(^±) & ve c -/+ :=] — oo, inf J_[ et J_ :=] sup J + , +oo[. Il 
est clair que I + H /_ =] sup J + , inf J_[. On demande que x± — 1 sur <^± tels 
que supp(x+) H supp(x_) = 0. Pour le comportement de x± à l'infini, nous 
imposons que x+ — 1 sur ] ~~ °°, m f ^+[ e t) X- — 1 sur ] SU P +°°[- Les x± 
s'annulent près de dl±. 

Considérons sur I + , une solution de l'équation P a e + = 0, 

% r ~ 

e+ := c+(a;h)exp(-- J g(y,a;h)dy). (2.8) 

Si nous choisissons le formalisme "BKW" , e+ s'écrit 

e + = [c + (a;h)e- ihf: + 9 * (y ' a;h)dy ^ f: + g{y ' a)dy , 

px 

=: a+(x, a; h)e^ +{x > a \ <p + (x, a) := - g(y, a)dy, 

J x+ 



avec une phase (p + G C°°, indépendante de h, qui vérifie l'équation eikonale 
d x ip + (x,a) +g(x,a) = 0, et une amplitude a + admettant un développement 
asymptotique en puissances de h, 

a + (x,a;h) ~ ^ a +ik (x, a)h k dans C£°(R), Va G fl 

Nous avons ip + (x + , a) = 0, ip' + (x + , a) = £ + G M, et 

Imv?'|(a; + (a),a) = -Im^(a; + , a) = > 0. 

La partie imaginaire de la phase est positive, cela restera vrai globalement 
sur I + : lmip + > 0. 

Grâce au choix de g vers — oo, pour C > assez grand nous avons 
||e + || i 2(]_ 0O) _ C ]) = 0(e~cK). H est alors possible, grâce au lemme de la phase 
stationnaire (que nous rappellerons après), de choisir c + (a; h) de la forme 

c + (a; h) ~ /T 1/4 (c° (a) + hc\(a) + ...)> 0, 
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avec 

a/4\ ï 



c o {a) = ^ -2lmg' x (x + (a),ay/^ 
{-lmg' x {x + ,a)) 1 / A 



ttV4 

tel que e+ soit normalisé dans L 2 sur /+. Ce qui implique, dans le formalisme 
"BKW" , que 

a + (x,a;h) = c + (a; K)e ih ^ 9 *^ dy 

= c+(a;h)(l-ih ! g°(y,a)dy + 0(h 2 )), 

J x+ 

\ïmg' x {x + ,a)\ 1 / A 



(tt/i) 1 ^ 



'x+ 

[1 + 0{h)) dans Cf 



Lemme 2.2 On désigne par e + G H SC {I + ) la solution normalisée dans L 2 (I + ) 
de (hD x + g(x, a;h))e + = sur I + . Nous avons la représentation asympto- 
tique suivante pour h <C 1, assez petit 

e+ ^ \Img' (x + g)|V* + 
(7m) 1 ' 4 

De manière analogue nous avons : 

Lemme 2.3 On désigne par e_ G H SC (I_) la solution normalisée dans L 2 {I ) 
de (hD x + g(x,a; h))*e_ = sur I_. e_ admet la représentation asymptotique 
suivante pour h <C 1, assez pefà£ 

e_ ~ (Im ^-^ ))1/ \ l + 0(fe))e-^-^r 

Nous rappelons ici le lemme de la phase stationnaire : 

Proposition 2.4 (Phase stationnaire) Soit a G C~(M). Supposons que 
G K — supp(a) et 

<p(0) = <p'(0) = 0, <p"(0)^0, Imy9>0. 

Supposons en plus que <p'(x) ^ sur K — {0}. Posons g{x) = (p(x) — 
{p"(0)x 2 /2. Alors, nous avons le développement asymptotique explicite, lorsque 
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h-ï 0, 

e i{p{x)/h a(x)dx ~ 

En particulier, nous voyons que le premier terme est (27ri/i) 1//2 y/'(0) _1 / 2 a(0). 
Pour une preuve on pourra consulter [19] p. 36. 

Proposition 2.5 Pour v G L 2 (I + ), v + G C, le problème de Cauchy (hD x + 
g(x, a; h))u = v, u(x + ) = 0, admet la solution unique u = Fv avec 

ii pu C 

\\F\\ L 2( I+ )-+H sc (I+) < -j=- 

Par ailleurs F admet un noyau intégral k vérifiant l'égalité 
\k(x, y)\ < 0{l/h) exp(-|x - y\/(CVh)). 
Preuve. F admet le noyau intégral 

k(x,y) = ^ e 4/;^;^ 1{i+ ^ i} _ '-e-yy^^h^y^y. (2.10) 

Nous montrerons après qu'un tel noyau vérifie l'inégalité suivante 

\k(x, y)\ < 0(l/h) exp(-|x - y\/(CVh)) (2.11) 

pour une constante C donnée. 

La norme L? de F est majorée (lemme de Schur) par 

/\ V 2 / r \ V 2 

|A;(x, j I sup / y)|rfar 1 . (2.12) 

Ce qui nous donne 

/C f C 

\k(x,y)\dy < —= et, sup / |fc(x, 2/)|cfe < — . (2.13) 

yh y ei+ J y/y 

Donc 

~ C 

\\F\\ L 2 {I+) ^ L 2 {I+) < — . 
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De la même façon on peut montrer que 

Il 2?|, C 

If f\\ L 2 {i+) ^ L 2 [i+) < —. 

Ce qui implique, en utilisant {hD x + g(x, a; h))Fw = w, que nous avons 

~ C 
\\F\\ L 2 (I+) ^ Hsc(I+) < — . 

Il nous reste à prouver l'estimation ( I2.1ip . Parce que x + est un zéro non- 
dégénéré de Img, et au comportement de Img vers — oo, nous avons pour h 
assez petit, et x G I+, (rappelons que g*(x, a; h) est uniformément borné) 

lmg(x, a; h) < — — (x — x + ) + C$h 2 si x > x + , (2-14) 
Co 

lmg(x, a; h) > — — (x — x+) — C\h 2 si x < x + . (2-15) 
Ci 

Ce qui implique, si nous considérons le cas x > x + , 



\k(x,y)\<le^y lm ^ h)d % x ^ 



h 

Pour (x — x+) > \/h, 12.161 s'écrit 



< L e -*&**-*+) - [y ~ x+) 'e^'l^^j, (2.16) 



\ -^h{{x-x+)-2{C hf)(x-y), ^ ^ 1 _ \ p -\x-y\/{CVh) 

et pour (x — x + 

Les autres inégalités nécessaires pour déduire (12.1 ip se calculent de la même 
manière. □ 

Proposition 2.6 Pour P a = hD x +g(x,a; h), v G L 2 (I + ), v G C le problème 

^ aU = V (2.17) 
R + u = v + 
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avec 



R + u : = (w,x+e+) = / u(x)x+(x)e + (x)dx, 
Ji+ 

admet une solution unique 

u = Fv + F + v + G H sc , 



ou 



avec 



1 1 

(e+,X+e+) D + 



D + = 1 + Q(e 



Ch 



et nous avons 



\F\\ L ^ Hsc = 0^), (2.18) 
\\F\\ c ^ Hsc = 0(1). (2.19) 



Preuve. Le problème 



P a u = 
R + u = v + 

admet la solution unique u = F + v + , tandis que le problème 

P a u = v 
R+u = 

admet la solution unique u — Fv :— (1 — F + R + )Fv. □ 

Nous désignons par L 2 comp (lJ) et H SCjComp (I_) respectivement l'ensemble 
des fonctions / appartenant à L 2 et à H SC {I_) respectivement telles que / 
est nulle dans un voisinage de inf I . 

Proposition 2.7 Pour P a = hD x + g(x, a;h),v<E L 2 comp (I_) le problème 

P a u + R-U- = v, 

avec 

:= -u_x_e_, -u_ G C, 
admet une solution unique dans H SCtComp x C, donnée par 

u = Gv, (2.20) 

u- = G_v := {x _ e \ e _ } (v,e-) =: ^-(v,e.), (2.21) 
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ou 

D_ = l + 0(e-à), 
avec, en tant qu'opérateur L^ omp (J_) — > H SCtComp (I ) : 

\\G\\l^h sc < 

Preuve. Pour x < rr_, nous avons l'unique solution de P a u = v, v G 
Ll omp (I-), qui est nulle près de inf J_ 

Ul (a;) = % - f X e-rJy~^ h ^v(y)dy := G x v(x), 

h J-oo 

alors que pour x > X-, nous avons 

u 2 (x) = % -\ X e-ï f y° {i ' a > h)d *v(y)dy := G 2 v(x), 

les deux exposants étant décroissants (à partie réelle strictement négative loin 
de X-) dans le domaine d'intégration. Afin d'obtenir une solution continue, 
il faut imposer 

• r+oo 

= Ul (a;_) - u 2 (x_) = -\ e-ï f y X ~ 9{ *' œ ' h)d *v(y)dy, 

h J-oo 

Donc, avec 

U- = G-V, 

nous pouvons prendre v — v — R-U-, car par construction 

(v — e_) = 0. 

Nous avons alors la solution 

u = Gv = G(I — R_G„)v, 
U- = G-V, 

où G est donné par G 12 dans les zones correspondantes. En observant que 
G a le noyau intégral 

k(x,y) = ^e-ïfy X °^ d % y < x < x _ } - j< rï,^^'^ Ur ^ 
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- qui est semblable au noyau intégral de F- et que lmg' x (x + ) ~ —ïmg' x (x_), 
nous pouvons utiliser les estimations du paragraphe précédent pour trouver 
que 



|G|U 2 (/_)->i 2 (/_) < 



C 

ainsi que 

C 

||G f || i 8(j_)_ >fl -„(j_) < -^=. 



□ 



Problème de Grushin global. Nous choisissons une partition de l'unité 
de M, ^-t G C°°(I±), ip+ + ip- = 1, telle que x± ~^ V'i, où ■?/>-< signifie 

supp O) n supp(l - 0) = 0. (2.22) 

Comme les x±i l es i } ± s'annulent près du bord de l'intervalle I±. On suit 
Hager |?| "Inverse global" : on commence par résoudre le problème sur I + 

(hD x + g(x, a; h))u = ip+v 
R + u = v + . 

Sur I + , on introduit 

ui ■■= (1 - X-)Fip+v + (1 - X-)F+v+. 

Donc, en utilisant x + -< (1 — X-)> e ^ i ) + ~^ (1 — X-) ( car X- ~^ ^-) 

R+ui = (m,x+e+) = v+, 

PaUx = 4> + v - [P a , X-)Fip + v - [P a , X-]F+v + , 

en rappelant que P a := hD x + g(x, a; h). 

Comme chez Hager [7j, on peut "corriger l'erreur" sur J_ en y résolvant 
le problème 

(hD x + g(x, a; h))u 2 + R-U- = ip-V + [P a , X-] F ^+ V + [ p a, X-] F + V +- 

Nous avons la solution 

u 2 = G(^v + [P a , X-}Fïp+v - [P a , X-}F+v + ) 
u_ = G-ty-V + [P a , X -]Fil>+v - [P a , X -]F + v + ). 

Et R + u 2 = car supp(u2) PI supp(x + ) = 0. 
Finalement nous obtenons : 
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Théorème 2.8 Pour tout a E il, 

v= (^hD x + g(x,a;h) ^ : ^ x c _> x c> 

est inversible d'inverse 

E E + 



E_ E. 



+ 



ou 



E = g(^_ + jx'-F^ + (1 - X-)Fi> + (2.23) 

E + = (l- X -)F + + G-x'_F + 

i 

E_ = + -X-F^ + ) 

Les normes vérifient 

\\E\\ L ^ Hîc = <D{±=), \\E + \\ = 0(1), \\E-\\ = 0(1), (2.24) 

\\E-+\\=o(Vhe-™y (2.25) 

Les opérateurs E, F + , G, G_ ont été définis au-dessus. De plus E + et E_ sont 
des opérateurs de rang 1, satisfaisant 

E + = ±^e + + 0(h~), (2.26) 

E- = (.,^e-) + 0(h~). (2.27) 

On peut déduire une formule exacte pour E_ + , 

}i C — i r x + 

E -+ = ~tT~tT~ ex P ( r / 9% a; h)dx). (2.28) 
%D+D_ h J x 

Finalement, en utilisant la relation (hD x +'g)~ l = E — E + (E V )~ 1 E_, on 

trouve l'estimation suivante pour la résolvante : 
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Proposition 2.9 On rappelle que 1/(9(1) > a < 0(1) appartient à un ou- 
vert Q, relativement compact et séparé de l'origine par une constante indépendante 
de h. Pour tout a dans Q, la résolvante vérifie pour h assez petit 

\\(hD x + g(x, a; h))^ ~ V + 0{h)) 

h 1 /' i {Lm g x {x + , a))* (1m g x {x_, a))* 

+ 0(4=), (2.29) 

où £ vérifie 

e :=-lm f g(y,a)dy. (2.30) 

Preuve. Il suffit de remarquer que E h s'écrit 

hh i 
\E-+(a)\ ~ —={lmg' x (x +1 a) Im^(x_, a)) 4 
\/7r 



1 r + 

exp(--Im y </(j/,a)dy)(l + 0(/O), 



et que 



□ 



L'opérateur d'Airy complexe sur KL Considérons l'opérateur A = + 
ix sur la droite réelle, muni du domaine D(A) = {u G H 2 (R)\ xu G L 2 (R)}. 
Nous pouvons montrer que pour chaque z G C, la résolvante satisfait (voir 
Helffer [II]) 

\\{A-z)- 1 \\ = \\{A-^z)- l \\. 



Proposition 2.10 (J. Martinet) Pour liez > assez grand, il existe deux 
constantes positives Cq et Ci telles que 

1 4 3 1 4 3 

C |Rez|"4 eX p -(Rez)2 < \\{A- z) _1 || < Ci|Re z\~ï exp -(Re z) 3 . 
o o 

nous voyons que 

J={D 2 X + ix- z)F~ l = i(-D x - ix 2 + iz). 
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Après conjugaison par l'opérateur unitaire u(x) i-» u(— x), nous notons (*4o + 
iz) := D x — ix 2 + iz ce nouvel opérateur. Nous avons ||(.Aq — z^) 1 1 1 = 
||(*4o — iRe2) _1 ||. Le changement de variable x = (Rez)^, nous permet 
d'identifier Aq + iz avec |Rez|Q, où Q est donné par (y a été remplacé par 

x) 

~-D~ — ix 2 + i. 

(Rez)t 

Nous nous retrouvons alors dans la situation précédente avec ïmg(x, a) < 0, 
plus explicitement Q est de la forme hD x + g avec g = g. En reprenant les 
notations précédentes, nous avons 

h = ~~~~~3i (2.31) 

(Re z) 2 

g{x,z) := -ix 2 + i, (2.32) 

x±(z) = ±1 vérifiant Im g(x±, z) = 0, (2.33) 

-lmg' x (x + (z),z) = 2x + (z) = 2, (2.34) 

+ Img' x (x-(z),z) = -2x-(z) = 2. (2.35) 

De plus 

1 f x ~ 1 Z" -1 2 
— — Im / g(x, z)dx — — — Im / g(a;, z)cfa; = — (— 2 H — )(Re^) 3 ^ 2 

h J x+ h Ji 3 

= ^(Re^) 3/2 , (2.36) 

et 

h 1/2 (lmg' x (x + ,z) ïmgf x (x-,z))* = h 1/2 (2x+(z)) 1/2 = 2 1/2 (Re z)~* . 
On déduit alors le corollaire 11.41 



2.3 Estimation de résolvante près du bord de l'image 
du symbole. 

On suppose cette fois ci que a > tend vers zéro pour le cas modèle 
à paramètre hD x + g(x, a; h) sur la ligne droite. On va se ramener au cas 
précédent par "scaling" . Rappelons que (lemme I2.ip 

x±(a) = W /2 (- ~ 2 ) ' + 0(a). (2.37) 
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Le changement de variable x = (a) 1 / 2 ?/, nous permet d'identifier notre opérateur 
avec aQ où Q est donné par 

h „ q{a x l 2 y,a) h 2 . , . 1/9 
cr/^a a cr 

~ + ^(a 1 V**) + ~ h 2 (a 3 g*(a l/2 y, a; /i)) où h = hjc? 12 . (2.38) 
a 

On a ainsi remplacé les points x + (a) et x_(a), par des points y+(a) = a~^x+ 
et y~(a) = a~ïx_, où y± x =)=1. Les points y± vérifient Im f(y±, a) = 0, 
où / est défini par f(x,a) := g ( a x ' a ) . Nous avons Im/^^O, 0) > 0, et 
Im/(x, a) x (x 2 — 1), Nous pouvons alors appliquer la proposition 12.91 avec 
h := h/ a 3 / 2 comme notre nouveau petit paramètre semiclassique. 

Proposition 2.11 Pour tout a, la résolvante vérifie pour h <C o?l 2 et a 
assez petit 

. yfn exp(ilm£(a)) 

ah^(hnf>(y + ,a))V±(hnf>( y „,a))V* { 

+ 0(^y=), (2.39) 
«V h 

où h = h/ a 3 / 2 , et i vérifie 



/ ~ —dy, 1mf(y ± ,a) = 0. (2.40) 

Jy- a 



Rappelons que f 

Ensuite, nous avons les égalités suivantes, 

1 «3/2 [*+ g(oï/*y, a ) a 3 / 2 g(y,a),_ 

dy = —r- / — Twdy 



et, 



h h J y _ a h J a i/2 y _ aa x l 2 

1 f x+ 1 
= J l j g(y,a)dy =: j-4>(a), 



«^(-Im /;(!/+, a)) 1 ^^^., a)) 1 



/4 



1/4.1/2 ( Im ^( al/ 2 y+, gO)j Im (g'x( al/2 y-> a )) \ 

a a 1 /» a i/8 

= /i 1/2 (Im^(x+,a))2(Im^(x_,a))3 

= fc 1/2 (^K,^}(p+))*(^fcM(p-))*- 

Nous pouvons ainsi donner une formule invariante pour la résolvante. 
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Corollaire 2.12 Soit g satisfaisant U.3\) et fli.^p . Pour tout a, la résolvante 
vérifie pour h <C a 3//2 avec a assez petit 



Il (MJ. + a; yggcpj|Im^ ^ + 0( - )} 

(2.41) 

où h = h /a 3 / 2 , et £o vérifie 

£ := Im / g(y,a)dy, lmg(x±,a) = 0. (2.42) 

Remarque 2.13 Puisque \a±(a)\ x as, nous avons les estimations sui- 
vantes : 



2 = « 3 / 2 



^t{Pq,P^}(p+) = -Im^(x+,a) x a 1/2 , 
^r{p^,Pa}(P-) = -Im^(ar_,a) x a 1/2 . 



\\{hD x + g{x,a-K)Y x \\ < -—^e^ 



Donc, sous la condition h <C a 3 / 2 <C 1 assez petit, nous avons 

C Co?lï 

Nous retrouvons alors la formule usuelle donnée dans [1] et [TB^ . 

Résumons le cas où a tend vers zéro, le théorème 12.81 traite en effet du 
point a > ^|yj : 

Théorème 2.14 Rappelons que (hD x — g~(x,a;h))~ 1 = E — E + Ez\E_, avec 
\E-\\ vWp(Ilm4(oQ) x h 

(2.43) 

/i, £ seront rappelés par la suite, et E = (9( ^-^ 1/4 ). i? + et i? + sont des 
opérateurs de rang 1, satisfaisant 

E + = 1 -^e + + 0(h™) 1 
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où D± = 1 + 0(h°°), et e+, e_ sont normalisés, et vérifient 
e + (x,a;h) = ( ^ + f f + ' (1 + OÇh))e^ ^ 

ou + = avec (7^{£ + g, £ + <?}(p+)) > 0. Nous avons une estimation 

similaire pour e_. Donc 



x (i + + e>(— 1_ ), 



où ft, := h/ a 3 ^ 2 , et £q (est une intégrale d'action) vérifie 

£ := Im / g(x,a)dx. (2.44) 

3 Cas général en dimension 1 

Soit p G S^IR^m), indépendant de h. On choisit un point zq appartenant 
au bord de l'image du symbole £(p) = p(T*M), qui est fermé g la 
condition d'ellipticité sur p et la croissance à l'infini de la fonction d'ordre 
m. Nous supposons que p est de type principal en z , c'est-à-dire 

p(p) -z = 0^ dp(p) ^ 0, 

et que p est d'ordre 2 en z . 

Pour simplifier, nous supposons que p _1 (z ) ne contient qu'un point po- 

En nous appuyant sur Hager ([QJ) : pour zq G <9£ nous avons p^ 1 (z ) = 
{po}, avec {p,p}(po) = 0. Soit Hi_^ p ^ le vecteur réel tangent à l'ensemble 

{p\ àip>P}(p) = °}- Donc 

(-Ti.To) 9 s ^ p(s) := exp( S if è{PïP - } )(p ) C {^fep} = 0} (3.1) 

est une courbe orientée, et ^{p,p} est positif à la gauche de cette courbe, 
et négatif à sa droite. Concernant la ligne fl3.ll) . nous avons évidemment 
p(0) = po, et le sous-ensemble de C, p ({|^{p, p} = 0}) , coïncide avec <9£ 
près de zq. 
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p(p(s)) est aussi une courbe orientée, puisque 

j- s (p(p(s))) = (H à{ ^ }P )(p(s)) = -{p, ±-.{p,p}}(p(s)) ? 0. 

De plus dlmp AdRep = ^■{p,p}dÇ i Adx, donc p — » p(p) préserve l'orientation 
à gauche de p(s), et la renverse à droite. Ce qui veut dire, qu'au voisinage de 
z , si p(s) est parcouru dans le sens positif, on voit que S se trouve à gauche 
de la courbe p(p(s)) G <9E. 

On note le chemin p(p(s)) par 7 :] — Ti,T [— > C. Un vecteur tangent à 
<9£ est donné par 7 (s), et nous savons que 

i(s) = -{p,^{p,p}}(p(s)). (3.2) 

Il existe un voisinage de zq où tout point z appartenant à S(p) s'écrit sous 
la forme 

z = 7 (s) + cm(s), < a < —, -T x < s < T , z = 7 (0), (3.3) 

pour une constante Cq > donnée, où n(s) est le vecteur unitaire normal 
à la courbe 7(5), orienté vers l'intérieur de n(s) vérifie les relations 

n(s) = ^]||fy|) ^ l n ( s )l — 1- P &r 1 & suite le vecteur unitaire tangent à 7(5) 

sera noté u(s), u(s) = j^jj et |u(s)| = 1. 

Puisque p est de type fini d'ordre 2 au point zq, nous pouvons alors 
supposer sans perte de généralité que p'^(po) 7^ et que p' x (po) = 0. Nous 
avons dans de nouvelles coordonnées (x, £) centrées en p(s) - changement de 
coordonnées symplectiques dépendant explicitement de s, conjugaison par un 
opérateur intégral de Fourier de plus le gradient de p en (0, 0) est un multiple 
d'un vecteur réel - 

P(P, s) - 7(s) = P(P, s) - p(p(s)) 

= u(s)( nK ' ^ + co(x, s)), avec co = 0(x 2 + £ 2 ), 

U{S) 

où lu est C°° en s. N'oublions pas que p va dépendre de s car les nouvelles 
coordonnées dépendent de s. Aussi p'ç(0, 0, s)/u(s) est réel. Puis, grâce au 
théorème de factorisation de Malgrange pour des fonctions C°°, il existe un 
ouvert V de R 2 x R + x] - 7\,T [ contenant le point (0,0,0,0) tel que dans 
V, nous avons la factorisation 

p(p, s)-z = p(p, s) - 7(s) - iau(s) 

pU0,0,s) 

= u ( s )( , \ — i + u(x, f , s) - ia) 
u(s) 

= u(s)q(x,Ç,a,s)(£ + g(x,a,s)), q(x, Ç, a, s) ^ 0, (3.4) 



23 



où g et q sont des fonctions C°°, g s'annulant au point (0, 0, 0). De plus pour 
x = £ = 0, a = et s non nul, nous avons p(0, s) — z = p(0) — j(s) = 0, 
impliquant g(0, 0, s) = 0. Par commodité, nous supprimons la variable s. 
Pour x = £ = 0, a — 0, nous avons alors 

<?(0,0) = 0, (3.5) 

g(0,0,0) = ^^, (3.6) 
u(s) 

9i(0,0) = 0. (3.8) 

Comme p' x (po) est nul, le rapport P ^° s '^ est réel et non nul, et sera noté a; 

a := = 9(°' ' )- Puisc l ue 

"4(0, 0) = <4(0, 0, 0) . + 2 ^(0, 0, 0) . + g(0, 0, 0)<4(0, 0), 
nous avons aussi, 

a^(0,0) =0/^(0,0). (3.9) 
Nous avons alors - rappelons que c'est pour a = - 
1 \u(s)\ 2 

= -aimera;, £), 

ce qui implique 

{p, ^{P,P}}(0, 0) = u(s){aÇ + u, -ïmaw' x }(0, 0) 
= - M ( S )a 2 Im^(0,0) 

= - M (s)a 2 aIm^(0,0). (3.10) 

Im<4(0,0)^0, (3.11) 

ï(s) 7(s) 3 |7(s)| 



On déduit alors que 

et que 

Im «4(0,0) 



u(s)a 2 a | 7 (s)| 3 ^(0,0) 3 

~ |feèfep}}(o,o)| 2 ^(0,0)3 



(3.12) 
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De l'équation (12. 3p qui donne le comportement de x±(a) dans le cas 
modèle hD x + g(x,a) avec l'hypothèse que g' a (0,0) = — i, nous obtenons 
qu'il existe un voisinage de z pour lequel nous avons pour tout complexe z 
appartenant à 

z£dZ(p), p- 1 (z) = {p+(z),p-(z)}, p:=(z,0, (3-13) 

où 

* ± M = t <o K (o, o, ,)|)V» ( |Img J M s)| ) V! + o(«), 

e est le signe de ^(0, 0, s)Im «^(O, 0, s). En remplaçant ^(0, 0, s) et Img^O, 0, s) 
par leur expression respective (13. 7p et (13. 12p . nous obtenons 

x ± (z) = XoT ea 1 / 2 \p'ç(p(s))\ ( 2 ) 7 +(5(a), (3.14) 

z = j(s) + an(s), a > 0, 
où e est le signe de {p, ^{p,p}}(p(s))/p'^(p(s)). Nous avons que 

±hp,p}(p±)>0. (3.15) 
1% 

On prolonge g sur M de la même manière que dans le cas modèle à pa- 
ramètre. 



On distingue, maintenant, deux cas : 

• 1er cas, p^ 1 (z ) ne contient qu'un point {po}> 

• 2ème cas, p^ 1 ^) ne contient que deux points {pi,p 2 }. 



3.1 Estimation de résolvante : cas 1 
Proposition 3.1 II existe deux fonctions 

q-J2^ khk e Sd(Vois(0,0),l), (3.16) 

fc>0 

9-Y,9kh\ 9k e C^MVois^O))), (3.17) 

telles que - nos coordonnées (x,Ç) sont centrées en p(s) - 

P(x, - z ~ q(x, £, a, s; h)#{£ + g(x, a, s; h))#q(x, £, a, s; h), (3.18) 

dans 5'(Vois(0, 0), m) uniformément en z (z dépend de s et a), où z appar- 
tient à un voisinage de z . Si l'on tient compte de la factorisation {p — z) = 
n(s)q(x,£)(£ + g(x)), nous pouvons choisir q égal à (n(s)g(x, £)) 1//2 impli- 
quant que gi = 0. 
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Pour la suite, par souci de clarté, nous supprimons la dépendance en a et s. 
Preuve. Avec p(p(z)) — z = 0, Pç(p(z)) ^ nous obtenons un voisinage V 
de p{zq) des fonctions lisses % et go définies respectivement sur V et 7r x (V) 
telles que 

(p-z) = ns)q(x,0) 1/ H^ + U^)ns)q(x^)) 1/2 

avec qo(x(zo),£(zo)) ^ 0, poO^o) = — ^(^o) On peut rajouter z aux variables 
et nous avons toujours une dépendance C°° de z, et les équations ci-dessus 
restent valables dans un voisinage de z . Nous regroupons par ordre de h les 
termes de la forme de composition asymptotique : 

=q^{x^)g N {x) + 2q (x^)q N (x 1 + 9o(x)) 
+ G N (q , . . .,q N -!,g , ■ ■ -9n-i, x, £, z). 

Etant donné que q est non nul, nous avons une équation de la forme 

où Gn ne dépend que deq , . . . , qN-i,9o, ■ ■ ■ 9n-i- H est possible de déterminer 
et qjsr inductivement (voir |8J), puisque le théorème de Malgrange nous 

prouve l'existence d'un voisinage de p(z ) et les fonctions g^ et q^ avec les 

propriétés demandées. En itérant, nous obtenons les séries formelles g et q. 
Pour le dernier point de la proposition, il suffit de remarquer que le terme 

en h d'une composition de Weyl symétrique s'annule. □ 

Nos variables p = (x,£) sont centrées en p(s) := (x(s), £(s)) ; si nous 
introduisons l'opérateur unitaire T s par 

T s u(y)-=e iy ^ )/h u(y~x(s)), 
nous trouvons que l'opérateur P := TgPT" 1 satisfait 

p = OPh(p) °ù p(p) = p(p - pO))- 

Il existe alors un symbole q pour les variables (x,Ç), a et s, vérifiant 

q = xf^Oït 1 )^) 5 )' 5 + ou X — Ca°(' K x,t(V)) est une troncature à 

support compact indépendant de z (donc de a et s) valant 1 dans un petit 
voisinage V de (0,0), tel que 

(q) w T s (P - z)T;\q) w = n(s)(hD x + g(x, a, s) + h 2 g*(x, a, s; h)) 

microlocalement près de (0,0), (3.19) 
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(remarque de J. Sjôstrand, pour y arriver, il semblerait qu'il faut aussi une 
rotation symplectique) uniformément pour z (qui est une fonction de a et s, 
voir (13. 3p ) dans un voisinage donné de z , 

g*(x, a, s; h) ~ g°(x, a, s) + hgl(x, a,s) + ... dans Cg°, 

uniformément en a et s. 

Par souci de clarté, la variable s est supprimée. On notera 

g(x, a; h) := g(x, a) + /i 2 g*(x, a; h). 

Nous pouvons trouver une partition de l'unité de M 2 , Xj £ S(R 2 , 1), j = 
0, 1 telle que (voir Faure-Sjostrand par exemple) 

Xj > 0, xl + Xi = 1, avec X o e C co (Vois(p )), Xi e C£°. 
Nous avons 

NI 2 = llXi^f + \\Xou\\ 2 ■ 
En utilisant l'ellipticité de (P — z)~ x on peut montrer que pour \a\ = 
0(/iln(i) 2 / 3 ) 

||(P - zy'vf < M 2 Jxov\\ 2 + 0(l)\\ X iv\\ 2 + 0{h)\\v\\\ (3.20) 
où M = ||(-P — ^) -1 Xo||- Ce qui veut dire que si M := max(M , 0(1)) alors 

iiop-^n < (M 2 +e>(/i))i 

Nous souhaitons maintenant connaître le comportement de Mo, 

M =\\(n(s))- 1 Tr 1 r(hD x + g(x,œ,h))- 1 rTsXo\\. 

Du théorème 12.141 et du fait que E (_ est scalaire et T s est unitaire, nous 

avons que 

1 v 



Mo = \EZH x \\q w E+E..q w \\ + 0( 



et nous obtenons une expression pour q^E+EJq" 1 , en fonction des quasi- 
modes : 

q w E + EJq w = ^L^(ç- e+ )(., ((ç)TV>-0 + 0{h°"). (3.21) 

Nous avons besoin d'un théorème de Melin, Sjôstrand sur l'action d'un opéra- 
teur pseudodifférentiel sur une fonction BKW avec phase complexe. Dès lors, 
puisque e+ et e_ sont des fonctions BKW normalisées, et microlocalisées près 
de p±, nous avons le résultat suivant : 
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Lemme 3.2 Nous avons la formule 

WçfE+E-çrW r r(l + 0(h)). (3.22) 

mp+)\ 2 \q(p-)\ 2 

Nous utilisons maintenant l'expression ( I2.43P de E h . Puisque (p — z) = 

n(s)q(£ + g), avec \n(s)\ = 1 nous obtenons que 

r— -Tri x (h^ + 9^ + 9}(p + )) ï = (^{p,P}(p+))K 
\ q ( p+ )\-2 2i 2i 

et similairement pour p_. Nous avons donc, pour tout h <C a 3 / 2 <C 1, assez 
petit 

M = vWjfa^ r( i + (h)) + 0(-=J—), (3.23) 

où /z. := h/ a 3 / 2 , et ^ (qui est une intégrale d'action) vérifie 

£o(a) := lui / lmg(x, a, s)dx (3.24) 



(remarquons l'ajout de s qui est important pour la suite). £o( a ) = £o(z) 
s'écrit aussi sous la forme 

rx + (z) 

lo(z) — Im / ip(x,z)dx, p(x, (p(x, z)) — z — 0, 

</lE_(z) 



Im / £cte, (7 relie p_ à p + ), (3.25) 

J-yCp- 1 {z) 



Im / ;cg?£. 



Pour que M ait que croissance tempérée en h 1 , il faut que 

(« 3/2 x) Im J X+ g(y,z)dy<Ch\n^, 
pour un C > quelconque. Soit si 



1\ 2 / 3 



a<0{l) ( Mn- ) . (3.26) 



On déduit alors le résultat suivant : 
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Théorème 3.3 Soit P un opérateur h-pseudodifférentiel de symbole indépen- 
dant de h. On suppose que (p, z ) vérifie les conditions M. 9) , M. 10]) , M.lty) . 
Un vecteur tangent à dT,(p) dans un voisinage de zo est donné par 7 (s). Donc 
tout point z de appartenant à un voisinage de zq peut s'écrire sous la 
forme 

z = 7 ( s ) + an(s), -Ta < s < T , a > 0, z = 7(0), (3.27) 

où n(s) = ry(s)/|7(s)|. Il existe un voisinage de Zç>, pour lequel tout point de 
z à l'intérieur de vérifie 

z i dE, p-\z) = {p + (z),p-(z)}, p± = (x±,e±), (3-28) 



avec 

±^{p,p}(p±) > 0. 

Il existe une constante T* (< Tq, T\) telle que pour toutes constantes Cq, Ci > 
0, il existe une constante C2 > telle que la résolvante (P — z) -1 est bien 
définie pour 



h 2/3 

et satisfait l'estimation 



7,2/3 1 1 

s\<T„ KaKC^hln-fl^ h < —, (3.29) 



\\(P - *ni ~ ^°P(fo°W) r x (1 + 0(h)) + 0(^J— 

où h = h/ a 3 / 2 , et £q (qui est une intégrale d'action) vérifie 



£o(z) := Im / Çdx, (7 relie p_ à p + ). (3.30) 

Dans le cas des opérateurs à coefficients analytiques, on peut pousser a 
jusqu'à une petite constante. 



3.2 Estimation de résolvante : cas 2 

Théorème 3.4 Soit P un opérateur h-pseudodifférentiel de symbole indépen- 
dant de h. On demande que pour tout z dans l'intérieur de S, p~ l {z) soit fini 
et que {p,p}(p) 7^ pour p G p _1 (z). 
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On suppose que p vérifie les conditions U.9\) . U.10\) . et que l'hypothèse 
U.lfy) soit remplacée par 

P~ X {Zq) ={P1,P2}, Pj = 

{p, {P,P}}(pi) ^ o, 

{p,{p,p}}(p 2 ) ^0. (3.31) 

sera l'union de deux domaines E 1; S 2 à bord C°° près de z dont les 
bords peuvent ne pas coïncider, £j = p(vois(pj)) . 

Soit Hxs Pt p} le vecteur réel tangent à l'ensemble {p\ ^{p,p}(p) = 0}. 
Donc 

(-Ï\,T ) 9s4 pj(s) := exp(sH_L {pp} )( Pj ) C {^fep} = 0} 

est une courbe orientée. Près de zq, pour tout z* G <9S(j»), nous avons 
p^ 1 (z) = {pi(s*), p 2 (U)}. p(pj(s)) est aussi une courbe orientée notée Jj(s) := 
p(pj(s)). Le vecteur tangent à <9£ dans un voisinage de z est donné par ji(s) 
ou par 72 (t). Il existe donc un voisinage de z tel que tout point z de se 
met sous la forme : 

z = Ji{s) + aini(s), —T x < s < T , a ± > 0, z = 71 (0), (3.32) 
z = 72(s) + a 2 n 2 {s) , -Ti < s < T , a 2 > 0, z = 7 2 (0), (3.33) 

n\ = «72 (s) /|7i(s) | et n 2 = i A f2(s)/\'ji(s)\. Alors pour tout point z de T, dans 
un voisinage de z , nous avons 

z i SE, p~\z) = {p l + (z),p l _(z),pl(z),p 2 _(z)}, 

où 

±^fep}(p±) > 0. 

Il existe une constante T* (< Tq,Ti) telle que pour toutes constantes 
Cq,Ci > il existe une constante C 2 > telle que la résolvante (P — z)^ 1 
est bien définie pour 

h 2/3 1 1 

|s| < T*, — — < a < Ci{h ln — ) 2 /3 , h < —, a = max(Q 1 ,a 2 ), (3.34) 
Co h C 2 

et satisfait l'estimation pour z G Si R S 2 : 

(P - Z)' 1 ~ sup - - V / " — -t X 1 + 0(h)) 

^=iAhy^{p,p}(p\))ï(^{p,p}(pi)) l J 
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pour z G Sx fi E2 



[p - *ni = £gg<%^M> r x (1 + (h)) 



ft 1/2 (èfeî'}(/'+)) i (èfep}(/'-)) J 



pour z G S' fl S 2 : 



(P - -exp(Ilm^)) x (1 + 

^ 1/2 (èfep}(p 2 + ))Mèferf(p-)) 3 

où h = h/ a 3 / 2 , et t (qui est une intégrale d'action) vérifie 

t (z):=\m ! £dx,{i relie p!_ à p\). (3.35) 

Preuve. Fort du cas 1, nous pouvons donner les grandes lignes de la preuve. 
Comme précédemment notre raisonnement donne donc deux opérateurs T s 
et T~ l tels que 

rT s (P - z)T; l r = ni (s)(hD x + g^x, a; h)) 

microlocalement près de (0,0), (3.36) 

gi = gi(x, a) + h 2 g*^(x, a; h), 

De même pour p 2 . Soient Vi, V2 des voisinages arbitrairement petits près de 
P\iPï- Il existe 0j G C^°{Vj) tels que si 

Po = P + Yl °i 

j 

alors PÔ 1 (zo) est vide. Soit pj = p + Ylk,k^j@k- Donc pj 1 (z ) = {pj}. Soit 
Pq = pQ(x,hD x ) et Pj = pJ(x,hD x ). P est uniformément elliptique. On 
suppose que = 1 près de po- Nous écrivons 

Ro(z) := J2( p j - *r% + (Po - - E^) 

j j 

Nous avons 

(p - z)Ro(z) =1 + J2( p - Pj)( p j - *r% 

+ (P + P )(Po-z)-\l-J20i), 



31 



où (P — Pj)(Pj — z) 1 9j = 0(h°°), (voir [8]) et aussi pour le dernier terme. 
Donc (P- z ) = Rq{z){1 + K)- 1 où K = 0(h°°). 

Nous pouvons trouver une partition de l'unité de M 2 , Xj G S(R 2 , 1), j = 
0,1,2 telle que 

Xj > 0, Xo + ■ ■ ■ + xl = 1, avec Xj -< Xo G C£°. 

Nous avons 

2 

INI 2 = XI llxHI 2 - 

j=0 

On peut alors montrer que 

2 

||(P - «J-^H 3 < J2 M !WxA 2 + O(l)\\xov\\ 2 + 0{h)\\vf, (3.37) 

3=1 

où Mj = \\(Pj — z)~ 1 Xj\\- Ce qui veut dire que si M := max(Mj, 0(1)) alors 

IKp-z)" 1 !! < {M 2 + o{h))ï. 

Nous souhaitons maintenant connaître le comportement de Mj qui est 
aussi donné par l'expression, 

M 2 =\\(n 2 (s))- 1 S; 1 q^(hD x + g 2 (x,a 2 ;h)y 1 q^S s X2\\, (3.38) 

S t a un comportement similaire à T s . 

De la section précédente, Mj, j = 1, 2 satisfait 



w v 7rexp(^Imf (z)) - 1 

^ 1/2 (è{^P}(P+)) 4 (àfo rf(P-)) 4 v « 1/4 

~ 3 /2 

où hj = h/aij , et £ J (qui est une intégrale d'action) vérifie 

f (z) := Im j £dx, (7 relie pi à p{). (3.40) 



□ 



4 Exemples 

Dans ce qui précède, le paramètre a représentait la distance du point z à 
la frontière du pseudospectre. Dans ce qui suit, pour les exemples, la frontière 
ne sera représentée que par des morceaux de l'axe réel ou imaginaire. Donc 
la partie réelle ou imaginaire de z sera notre distance à la frontière, soit notre 
paramètre a. On préfère ici changer et utiliser la lettre y à la place de a. 
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4.1 Cas 1 : L'opérateur cubique non- autoadjoint 

On considère l'opérateur non- autoadjoint sur la droite réelle 

V = (D x ) 2 + ix 3 . 

L'image du symbole principal est le demi-plan Rez > 0. Le spectre est 
composé de valeurs propres discrètes situées sur l'axe réel positif (voir Tre- 
pheten, Embree [E]). Le pseudospectre est symétrique par rapport à l'axe 
réel. Pour étudier les courbes de niveaux de la résolvante, il suffit donc de 
considérer celles situées au-dessus de l'axe réel, p désigne le symbole de V . 

Pour z au dessus de l'axe réel et appartenant à <9£(p) = iR, nous avons 
p- 1 (z) = {p} = {(Imz)3,0} et 

{p,^fep}}(p) = -2x3 2 (|Im^|) 4 /V0. 

Nous sommes donc dans le cas 1. 

Pour tout z appartenant à l'intérieur de nous avons 

où 

P±0) = ((Imz)5,iF(Rez)5), 

et 

^{P,P}(P±) = (lmpçPx)(p±) = -2 x 3£x 2 . 

Le changement de variable x = (Imz)^x, nous permet d'identifier (V — z) à 
(ïmz)V où V s'écrit 

/ f „ \i , . 3 Rez . , 1 Re^ . . 
( -U x ) + ix z, avec h := g-, y := . (4.1) 

(Imz)« lmz (Imz)ê imz 

h est notre paramètre semiclassique et y est choisi petit. On choisit y > 0, 
puisque l'on se restreint aux courbes situées au dessus de l'axe réel. Soit </?(£) 
la fonction qui vérifie 

e + ^<p(O 3 -y-l = 0■ 

Nous avons 



r-vv r-w 1 

/ <pdÇ= (l-iy + i£ 2 )3de 
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Puisque y est petit, un développement limité donne 

Im / (l-iy+ ie)îdÇ = -ys + £% 5 / 2 ). 
■/^ y 

Ensuite, remarquons que 

= (-3x 2 e-3x 2 0(l,-2/ 1/2 ) = 6v^, 

il en découle que 

Tr-iPiP} (p-(y +i )) = Q Vv- 



Proposition 4.1 Pour Imz > 1 et y := Sa* ) : = — L_ 5 _ ) ^ e / s 



Im2 ' (Im2)« 



1 / ) 

3 



h « y 3/2 « 0(/iln-) sorf, (Im2:) 4/9 < Rez < Imz, 

nous avons de manière précise, 

v^rexp ( ±Im f~/"(l - iy + ix 2 )^dx) i 

11(7? _ z }-i|| ^ \1 il L + 0( - 

" l ^ " 6V2(W)i/3(Re*)V4 l (Re2) i (Imz) 

on frzen, sous une forme simplifiée 

./?cxpr 4 (Re2)3/2 i cr (Rez)5/2 n 1 

,Wll V Pl 9(Im^) 2 / 3 + U l(Im Z ) 5 ^ 1 , 

11 6V 2 (I m ,)V3(Re,)V4 + 0( (Rez) | ( i mz) | ) - (42) 

Lignes de niveaux de la résolvante. On se sert de l'expression (14.2p . 

1 



\\(v- Z )-'\\ = <m 



On obtient alors 



_ 4 (Rez) 3 / 2 (Rez) 5 / 2 / y^F 1 \ 

116 ~ 9(Imz) 2 / 3 + ( (Imz) 5 /3^ + 11 (R ez )i/4(l m 2 )V3 J • t-J 
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(14. 3 p se met sous la forme 

= (A)2/3( Imz) 4/9 (l n ( Re2 (I mz )4/3^ 2 / 3 + 0(^2)4/9) 

,(Re^ 5 / 3 



Donc 



Re2 =(A) 2 /3(i mz )4/9 rin(Imz) 16 / 9 (ln((Rez)(Imz) 4 / 3 ) 2 / 3 ) 2/3 
16 

+ 0((lmz) 4/9 ). (4.5) 
Puisque Rez < Imz et que ln(t ln£) = (1 + o(l)) lnt, on obtient 

Rez = (1 + o(l))(^^) 2 / 3 (Im2) 4 / 9 (lnImz) 2 / 3 . (4.6) 
16 ■ 9 

Résumons : 

Proposition 4.2 La partie des lignes de niveaux de l'opérateur non- autoadjoint 
V se situant en dessous de l'axe réel, admet la représentation asymptotique 
suivante lorsque Im z — > 00 

Rez = (1 + o(l))(Imz) 4/9 (lnImz) 2/3 . 

De manière plus précise, si l'on considère la ligne de niveaux - = \\ (V — z)^ 1 \\ 
alors nous avons 

Rez = (l + o(l))(^) 2 / 3 (Imz) 4 / 9 (ln( ^ mZ ^ )j * . 

4.2 Cas 2 : Oscillateur harmonique non-autoadjoint 

Maintenant, considérons l'opérateur harmonique non-autoadjoint, 

Q = {D x ) 2 + ix 2 (4.7) 

sur la droite réelle. 

E(p) est le premier quadrant supérieur, et <9E(p) = R + U iR + . 

Les lignes de niveaux de la résolvante sont symétriques par rapport à 
l'axe e î7r,/4 M. Intéressons-nous ici au comportement asymptotique des lignes 
de niveaux. Au regard de la symétrie axiale, il suffit de considérer la partie 
des lignes qui se situe en dessous de l'axe de symétrie. Nous avons le résultat 
suivant de K. Pravda Starov : 
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Proposition 4.3 Si est la symétrie axiale par rapport à l'axe e ln > R. Il existe 
C > tel que nous avons 

{z e C; Rez > u et < Imz < C{Rez) 1/3 - e}U 

Si ({z G C; Rez>u etO<lmz< C(Re z) 1/3 - e}) C cr t (Q) c , 

a e désigne le e -pseudo spectre de Q. 

Pour tout z dans l'image du symbole g, nous avons 

q -\z) = {pï_( Z ),p> + (z)\j = l,2}, (4.8) 
où^{g,ç}(p ± )>0. (4.9) 
et si z un point du bord, appartenant à R + alors 

g- 1 (^o) = {(eo,0)}n{(-eo,0)}, £o>0. 
Nous nous trouvons donc dans le cas 2. Dans un voisinage de z nous avons 
pl(z) = ((Imz)3,-(Itez)5), pi(^) = (-(Im^,-(Re^), 
P+(<z) = (-(Im^, (Rez) 3), p 2 _(z) = ((Imz)5, (Rez) 3). 

Le changement de variable x = (Rez)2x, nous permet d'identifier Q — z à 
(Rez)Q où Q s'écrit 

1 ^ V • 2 -, .Imz , 1 Imz , rt N 

.Dj; + ix - 1 - z— — , avec h := - — , y := - — . (4.10) 



Rez J Rez Rez Rez 

h est notre paramètre semiclassique et y est choisi petit. 
Nous avons 

£ 2 + ix 2 - 1 - iy =(£ - yjl + iy - ix 2 )(£ + yjl + iy - ix 2 ), 

=:(ç-<Pi(x,y))(t-<P2(x,y))- (4- 11 ) 

Nous obtenons 

rVv r-Vv rVv . 

/ ip 1 dx = ip 2 dx — I yjl + iy — ix 2 dx. (4-12) 

J-Vv J \/y J-y/v 

Nous souhaitons calculer cette inté grale. 

Soit F(x) = | arcsinx + x ^ 1 ~ x ' 2 ; = y/1 — x 2 , nous avons 



cfe \y/TTvy) ~ y ~ l + iy VTTJy' 
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Donc 

d „ / e'ïi 



e 4 (1 + iy)—F I , I = y 1 + iy — ix l . 



dx \ y/1 + iy 

Si bien que 

-y/y 



i7T_ , c 4 -\/y 

e 4 (1 + iy) arcsinf— v ), (4-14) 
V 1 + «y 



un développement limité donne 



V / ï+^v / ^+^ / 1 . i/ 3/2 + (% 5/2 ). (4.16) 
6 (l + iy)3 



Remarque 4.4 On aurait pu aussi traiter l'intégrale en remarquant 

que 



V x 2 + a = — (x v x 2 + a + a ln(x + \/x 2 + a)), 
ce donne 

, i y^+v^ 

/ ^ 2 = -e 4 (î _ y ) ln -. 

JVy 2 "v^ + vî 

(9n retrouve M -HP V ar ^ a formule : 

arcsin(z) = —i ln (iz + (1 — z 2 ) 5 J . 
Finalement, en prenant la partie imaginaire dans (14. 16p . on obtient 

Im <pi dx = Im ip 2 dx = V 2 + V 2 = \y^ + 0{y^) 

J-Vv Jy/y ° 6 

Ensuite, remarquons 

^{^, q]{p\{l + iy)) = ^{g, q}(y/y, -1) 

= (-2^-2^)(-l, v ^)=4 v ^, 

^{ç, + iy)) = 4Vy- 

On déduit alors le résultat suivant : 
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Proposition 4.5 Pour Rez 1 et y := h := 

/i < y 3/2 < 1 sozi (Rez) 1/3 < Imz < Rez, 
nous avons de manière précise, 

^ V^rexp^Im(e-T(l+i î/ )axcsin(^5))J 



\\{Q-z) 



2(Rez)V4(i mz )i/4 

+o{- 1 



(Imz) * (Re z) * 
on bien dans une forme symplifiée, 

./^r::nr 2 (Im2)3/2 I Q( {Iraz)5/2 )) 1 

ii (g - ,m ~ v P 3(Re : ) ;;: :r )3/ + o( 



2(Re^) 1 /4(i m2 )i/4 v (Imz) 1 /4(Rez) 1 /* y- 

(4.17) 

Lignes de niveaux de la résolvante. On se sert de l'expression (14.17p . 

1 



On obtient alors 



_ 2 (Imz) 3 / 2 (Imz) 5 / 2 /y^F 1 \ 

" 3 (Re Z )V2 + ° ( ( Re z f/2 ' + ln \~ (Re z)V4(i m z )i/4 J ' ^- 18 > 

En oubliant le facteur logarithmique dans ( 14.18p . on aurait 

2(Imz) 3 / 2 ^,(Imz) 5 / 2 N 2(Imz) 3 / 2 , ^lmz„ , 



Soit 



ainsi 

Imz= (Rez(-lne) 2 ^) 1/3 (l + 0(j^)). (4.20) 
(I4.18P se met sous la forme 

A 2 /3m m/s A /2(Re^) 1 / 4 (Im^) 1 /4 \ 2/3 (Imz) 5 /\ 

= (^(Re*) 1 /» (ln(Re.Im,)) 2 / 3 + ©((Hez) 1 /») + 

(4.21) 
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Donc 

Imz =(-) 2/3 (Re^) 1 / 3 (ln(Rez) 4 / 3 (ln(RezIm^) 2 / 3 ) 2/3 
6 

+ C((Re^) 1/3 ). (4.22) 
Puisque Imz < Rez et que ln(tlnt) = (1 + o(l)) lnt on obtient 

Imz = (1 + o(l))(— ) 2/3 (Rez) 1/3 (lnRez) 2/3 . (4.23) 
6 ■ 3 

Résumons : 

Proposition 4.6 La partie des lignes de niveaux de l'oscillateur harmo- 
nique non- autoadjoint se situant en dessous de la droite e î7f//2 M admet la 
représentation asymptotique suivante lorsque Re z — > oo 

Imz = (1 + o(l))(— ^(Rez^lnRez) 2 / 3 . 
6 ■ 3 

De manière plus précise si l'on considère la ligne de niveau - = \\(Q — z)~ 1 \\ 
alors nous avons 

Imz = (1 + o(l))(^) 2 / 3 (Res) 1/3 (in ^f^ ^ ^ 

On remarquera que les lignes de niveaux de la résolvante ne se rapprochent 
pas à l'infini puisque partant de (14.21p . nous avons 

Imz = (-) 2/3 (Rez) 1/3 (ln(Re z) 4 / 3 (ln(Re zlm z) 2/3 - lne) 2/3 

(Imz) 5 / 3 
+ { (Rez)^ ), 

. Co(Re2 )V3 (lnRe , )2 /3( 1 + Cij z^i +0( jH_ )2 ) 

(Imz) 5 / 3 
+ l (Rez) 2 / 3j ' 

=C Q (Rez)^(lnRez)^ + C' {Rezf/ 3 ~\ /3 +... 

{mRez) 1 ' 

Le second membre tend vers l'infini pour un e donné. 

Projection spectrale. Ici nous nous intéressons au projection spectrale de 
l'oscillateur harmonique non autoadjoint. 
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Théorème 4.7 (Davies, Kuijlaars) Soit l'opérateur H = D 2 + z A x 2 pour 
un complexe z = e ld . Les valeurs propres sont notées \ n = z 2 {2n+\). On note 
N n>z la norme du projecteur spectral ^ j D ^ x ^{z — H)~ 1 dz. Si < 6 < 7r/4 
alors 

lim n~ l ln N njZ = 2Re {f{r{6)e ie )) 



ou 



f(z) = zVz 2 - 1 + \n(z + Vz 2 - 1), 

et 

r{6) = (2cos(2^))- 1 / 2 . 

Nous allons redémontrer ce théorème avec nos estimations. 
Preuve. Le changement de variable x = (2n + \) l l 2 x permet d'identifier 
H - \ n à (2n + \){{hD x f + z 4 x 2 - z 2 ), avec h = (2n + Compte tenu 
de qui précède, nous avons 



lim \nN nz = -hm [ + ^(z 2 - z 4 x 2 ) dx = -hm [ + ^(1 - z 2 x 2 ) 

n^oo h J x _ h J x _ 

= ^Re f + ^/(z 2 x 2 - 1) dx, 

h Jx- 



dx 



ou 



Imz 4 y sin(4#) 

= (2cos(20)) 



sin(2#) \ _ m »»„mn\\-i/2 



i 



2sin(2fl) cos(2fl) 
x _( z ) = -(2cos(2^))- 1 / 2 . 

Etant donné la remarque 4.4, il est facile de conclure. □ 



4.3 Opérateur d'advection-diffusion 

Ici nous calculons la résolvante d'un opérateur d'advection-diffusion au 
point 1 + 

On considère l'opérateur non-autoadjoint sur le cercle 

L = — sm(x)(hD x ) 2 — ihD x . 
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Fort des informations de [TB], nous savons que l'image du symbole à l'infini 
est comprise entre les deux paraboles liez = ±(Imz) 2 . Nous allons montrer 
que c'est exactement la frontière de son image. Pour cela, remarquons que 

t x {z) = {(arcsin(--5^-), -lmz)} U {(n - arcsin(- R<32 ' ), -Imz)} 
[lmz)' s {lmz) z 

= { P+ }U{p_} (4.24) 
Nous avons 

^{£,£}(p)=ecos(x). (4.25) 

Par souci de simplification, nous notons 7 := n™^a ■ Puisque cos(arcsinx) = 
vl — x 2 . nous voyons que 

±{l,i}(p + ) = Qmz)*y/ï=ï*>Q, 
±:{t,l}(p-) = -(Imz) 2 v^7 < 0. 

En utilisant l'annulation du crochet de Poisson sur le bord de l'image de l 
nous obtenons : 

Lemme 4.8 Soit L = — sm(x)(hD x ) 2 — ihD x de symbole £ alors£(T*S 1 ) est 
bornée par les paraboles 

liez = -(Imz) 2 et liez = +(Ini2) 2 . 

De plus, nous affirmons que les points du bord de £(£) sont tous d'ordre 2, 
à l'exception de l'origine. 

Preuve. Nous allons calculer le deuxième crochet pour vérifier que les point 
du bord sont d'ordre 2 

R ÔXp) ={- sin(a;)£ 2 - i£, £ 2 oos(ar))} 

=2£ 3 (sinx) 2 + 2£ 3 (cosx) 2 

+ iÇ 2 sin x 
=2£ 3 + ^ 2 (sinx). 

Dans le cas où Imz ne s'annule pas, il est clair que le bord est d'ordre 2. □ 

Soit le point (1+z)— ae l6 6 g] — 7r/2, 7r/2[ où cos 
Nous partons pour a petit dans la direction donnée par la normale à la courbe 
dH{£) au point 1 + z. 
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Nous avons les formules exactes plus les développements limités suivants 
(en utilisant arcsin(l — x 2 ) ~ | + \/2|a;| + C(|x| 3 )) : 

. 1 — a cos 9 

x+(a) = arcsm(-- r-^ô), 

(1 — asm9)' s 

= -^ + a 1/2 (V2x (-cos9 + 2sm9)) 1/2 + 0(a) 



7Ï 



^ + a 1 / 2 ( v / 2^5) 1/2 + 0( 

, , . . 1 — a cos 9 

x_{a) = 7i — arcsm( 



a 



1 — a sin 9) 2 



En utilisant la formule 13.141 on aurait retrouvé 
Pour £± nous avons 

£± ( a ) = — 1 + a sin 
Avec les données exactes de x± et £± nous retrouvons 

i{a}( P± )xa^. 

Pour l'intégrale d'action, la fonction implicite (p(x, z), vérifiant tp(— 7r/2, 1+ 
i) = — 1, de l'équation —(sin :r)£ 2 — z£ — z = est 

1 x -* + (-l-^sin(x))V 2 ; 

2 sinx 

et 

— z cos a; 



2(sinx) 2 

2z(cos x)(— 1 — 4zsinx)~ 1 / 2 — (cosx)(— 1 — 4zsin:r) 1 / 2 



+ 



(sin x) 2 



Comme ip' x (— tt/2, z) = 0, et <^(— 7r/2, z) ^ on peut retrouver 

r x + 

£ {z) = ïm <p{x,z)dx x a 3/2 . (4.27) 



Nous avons alors redémontré dans le cas de l'opérateur d'advection-diffusion, 
la formule que nous connaissions [16] : 



42 



Proposition 4.9 Soit 9 tel que cos 9 = —1/y/E, sin# = 2/y/E. Pour h <C 
a 3/2 < i nous avons 

1 expfiflm [<p(x, (1 + i) - ae id ))] x +) h 
+ C?( 



on 



avec 



1 -i + (-1 - 4zsin(x)) 1/2 

2 sinx 



. 1 — a cos 

x+(a) = arcsm(-- r-^ô), 

(1 — asmé*) 2 

1 — et cos 

x_(a) = 7T - arcsm(-- r-^â). 

(1 — asmO) 2 



Ce qui donne, après simplification 
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